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The collective behavior of spins in a dilute magnetic semiconductor is determined by their mutual 
interactions and influenced by the underlying crystal structure. Hence, we begin with the atomic 
quantum-mechanical description of this system using the proposed variational-perturbation calculus, 
and then turn to the emerging macroscopic picture employing phenomenological constants. Within 
this framework we study spin waves and exchange stiffness in the p-d Zener model of (Ga,Mn)As, 
its thin layers and bulk crystals described by the spds* tight-binding approximation. Analyzing 
the anisotropic part of exchange, we find that the Dzyaloshinskii-Moriya interaction may lead to 
a cycloidal spin arrangement and uniaxial in-plane anisotropy of diagonal directions in thin layers, 
resulting in a surface-like anisotropy in thicker films. We also derive and discuss the spin-wave 
contribution to magnetization and Curie temperature. Our theory reconstructs the values of stiffness 
determined from the temperature dependence of magnetization, but reproduces only partly those 
obtained from analyzing precession modes in (Ga,Mn)As thin films. 



I. INTRODUCTION 

Dilute magnetic semiconductors, such as (Ga,Mn)As, 
form a class of materials with many outstanding 
propertiesii^ and functionalities^^, A number of these re- 
sult from the complex structure of the valence band, 
which hosts holes mediating magnetic order between 
spins localized on transition metal impurities. The 
strong force which lies behind this mechanism is the 
carrier-mediated exchange interaction. It gives the low- 
lying energy states of such a system the character of 
spin waves, which contribute to both its equilibrium 
(e.g. spontaneous magnetization and Curie temperature) 
and nonequilibrium magnetic properties (ferromagnetic 
resonance and relaxation))^ At the same time, it allows us 
to replace the atomic quantum-mechanical description of 
the system with the classical, continuous micromagnetic 
theory, where it appears in form of exchange stiffnessi^ 

Quite generally, micromagnetic properties of any fer- 
romagnet are determined by the magnitudes of the 
exchange stiffness and magnetocrystalline anisotropy?^ 
The first describes the exchange energy associated with 
nonuniform distributions of local directions of magneti- 
zation. The other is the energy needed to change the 
total magnetization direction with respect to the crystal 
axes, which involves the competing crystal-field and spin- 
orbit interactions. While the main part of the exchange 
energy is isotropic, a consequence of its electrostatic ori- 
gin, the relativistic spin-orbit coupling can create its 
small anisotropy, namely the dependence on the crys- 
talline orientation of magnetization. In zinc-blende bulk 
crystals and thin layers with broken space inversion sym- 
metries, it leads to many interesting effects, related to 
the anisotropic^ and Dzyaloshinskii-Moriyai^i^ exchange, 
which have not been hitherto studied in (Ga,Mn)As. 
Since the total spin is no longer conserved, they are in- 



dispensable when considering any usage of spins for in- 
formation storage and processing. 

The spin-wave spectrum and the isotropic exchange 
stiffness in bulk (Ga,Mn)As were computed by Konig, 
Jungwirth and MacDonaldfi^ and by Brey and Gomez- 
Santosii within the p-d Zener model employing the six- 
band k ■ p Hamiltoniaufi^r— which neglects the inver- 
sion asymmetry specific to the zinc-blende lattice. It was 
found that the actual magnitude of the exchange stiff- 
ness is much greater when one takes into account the 
complex structure of the valence band, as compared to 
the case of a simple parabolic bandJ^i^ This, as well 
as the highly anisotropic Fermi surface, were shown to 
explainiiii why the mean-field approximationi^ is so ac- 
curate in (Ga,Mn)As. It was also found that the val- 
ues of exchange and anisotropy energies obtained within 
the same formalism describe quantitativelyi^ the width 
of stripe domains in films with perpendicular magnetic 
anisotropy. 

More recently, Bouzerar— employed a self-consistent 
local random-phase approximation in order to evalu- 
ate the spectral density of spin-wave excitations in 
Gai_2;Mna;As. The magnitudes of spin-spin exchange 
integrals J{r) were obtained from first principles com- 
putations within the local spin-density approximation 
(LSDA) and tight-binding linear muffin-tin orbital ap- 
proach neglecting the spin-orbit interaction. The the- 
ory allows to treat disorder and thermal fluctuations, 
and shows that the calculated spectral density has well- 
defined maxima up to about one half of the relevant De- 
bye wave vector qd = (24a;/7r)^/^7r/a, where a is the lat- 
tice constant. This made it possible to determine the 
spin-wave dispersion uj{q) in the range < g < 90/2, 
from which the magnitude of spin-wave stiffness was 
obtainedfi^ 

Experimentally, Potashnik et alr^ analyzed the tem- 
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perature dependence of magnetization in a series of 
Gai-ajMn^^As samples, which provided the values of spin- 
wave stiffness from the T'^/'^ Bloch law. In later ex- 
periments, the stiffness was determined by examining 
spin precession modes excited by optical pulses^ and 
under ferromagnetic resonance conditions.—"— The val- 
ues obtained for some films with thickness greater than 
120 nm, either as-grow n^^'^^ or annealed)^ are in good 
agreement with those predicted by BouzerarJ^ How- 
ever, the values for thinner film o^^'^^ or another series 
of annealed sample s were about three times smaller. 
Similarly small magnitudes of spin-wave stiffness were 
found by analyzing the domain structure of annealed 
Gao.gaMno.orAsj^ The experimental workai^r^i demon- 
strate that spectral positions of spin- wave resonances are 
strongly affected by the character of spin pinning at the 
sample borders. The same effect is caused by magnetic 
anisotropy changes along the growth direction, particu- 
larly strong at the film interface and surfaccj^ Therefore, 
theoretical predictions concerning spin-wave excitations 
in both thick and thin layers may provide a useful guide 
to better understanding of magnetization dynamics in 
real samples. 

In this paper we investigate spin waves and related 
micromagnetic constants in ferromagnetic (Ga,Mn)As, 
as described by the p-d Zener modelJ^ The validity of 
this model is supported by photoemission experiments^ 
and ah initio computations, in which the inaccuracies of 
the LSDA are partly reduced by self-interaction.— Fur- 
thermore, the magnitude of low-temperature quantum 
corrections to conductivity indicates that in the con- 
centration range relevant to ferromagnetism, the den- 
sity of states assumes values expected for valence band 
holeSi ^^'^^ The carrier band is formed from p-type-like 
states of the GaAs semiconductor structure, which are 
mostly built from the anion 4p orbitals, but has a rather 
large d component as the result of a strong sp-d inter- 
action. We describe it by the spds* tight-binding ap- 
proximation for thin layers and, by applying periodic 
boundary conditions, for bulk crystals^i""— Our analy- 
sis of the anomalous Hall effect in (Ga,Mn)As in Refs.lSlI 
and [s^ provides a thorough comparison of this method 
with other band structure models<^ One of its advan- 
tages is that it captures the inversion asymmetry of 
the zinc-blende lattice, which produces the Dresselhaus 
spin splitting of the conduction band)^ and addition- 
ally the structure inversion asymmetry in thin layers, 
which creates the Bychkov-Rashba spin splitting)^ The 
growth-induced biaxial strain is included by changing the 
atoms' arrangement, according to the strain tensor val- 
ues: Exx = £yy ~ 5a/ a and Ezz ~ —'^^12/ chExx, where 
5a is the strain-induced change of the lattice constant, 
and C12/C11 = 0.453 is the ratio of elastic moduli. Also, 
the on-site energies of the d orbitals depend linearly on 
the strain tensor values.— The sp-d exchange coupling 
is modelled using the virtual crystal and mean-field ap- 
proximations, while taking into account the appropriate 
weights of Ga and As orbitals in the wavefunctions close 



to the center of the Brillouin zonei^i^ In this way, spin 
polarization of individual Mn moments is replaced by a 
molecular field, which creates a fc-dependent Zeeman-like 
splitting of the host bands. For heavy holes in the F point 
it is given by A = xhqSP, where x is the fraction of cation 
site density uq substituted by Mn forming a spin S ~ 5/2, 
and /3 = — 54 meVnm^ is the p-d exchange integral 
The validity of this approach can be questioned if the 
magnetic ions produce bound states^^ but the metallic 
character of the carrier states of interest here means that 
this does not occur due to many-body screening. At the 
same time we neglect the effect of this screening on po- 
tentials produced by magnetization fluctuations, which 
is justified as long as A is smaller than the Fermi energy, 
so that spin and charge density fluctuations are decou- 
pled. A part of Mn atoms form unintentional defects 
such as intcrstitials (which can passivate single substitu- 
tional spins) and antisites. Both being double donors, 
they signiflcantly lower the hole density, which can be 
partly remedied by removing the intcrstitials by post- 
growth annealing. 

The paper is arranged as follows. In Sec. [H] we use the 
proposed variational-perturbation calculus to describe 
the system of lattice spins interacting via hole carriers, 
and its spin-wave excitations. In the micromagnetic the- 
ory, the system is described by the set of constants related 
to magnetocrystalline anisotropy and exchange energy. 
Section III CI provides quantitative results on the spin- 
wave stiffness expressed as the dimensionless parameter 
Dnor, comparing it between the k ■ p and spds* tight- 
binding models. In Sec. IIIDI we derive the spin-wave 
contribution to magnetization and Curie temperature. 
Section [HT] compares our theory to related experimen- 
tal findingsJi"— The last section of the paper contains 
the summary of our work. 



II. THEORETICAL MODEL 

In this section, we demonstrate a physically transpar- 
ent perturbation- variational method of treating the sys- 
tems in question. We use it to find the effective Hamil- 
tonian of lattice ions interacting through hole carriers, 
written in the interaction representation of creation and 
annihilation operators. Next, we calculate the disper- 
sion dependence of its low lying energy states by trans- 
forming the Hamiltonian to the spin-wave representation. 
These two pictures are associated with different sets of 
phenomenological parameters describing the macroscopic 
system, whose properties we shall investigate later on. 



A. Microscopic picture 

We consider the ferromagnetic phase of a system con- 
sisting of P carriers and N magnetic lattice ions, de- 
scribed by the Hamiltonian Ho a-nd coupled by the p-d 
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exchange interaction %' , 



N 



H = Ho+H' = Ho + ^^/?/(r,-R,)s, -Sj , (1) 

i=i i=i 

where and Sj are the i-th carrier's and j-th ion's 
spin operators, while and are their respective po- 
sitions. The strength of the p-d exchange interaction be- 
tween these two spins is described by a smooth function 
jSIiri — Rj), locahzcd around the j-th magnetic ion. In 
the absence of external fields, "Ho depends on the carriers' 
degrees of freedom only.— 

The dynamics of magnetic ions coupled to the sys- 
tem of hole carriers requires a self-consistent description, 
which takes into account how the holes react to the ions' 
magnetization changes. Therefore, we use the Lowdin 
perturbation method specifically adapted for multiparti- 
cle HamiltonianS)^"— to derive an effective Hamiltonian 
^cft £qj, ^pj^g only. 

We choose the multiparticle basis states of % as M ®T. 
The ion part M is an eigenstate of the unperturbed ions- 
holes system ([1]), and the hole part F is a Slater determi- 
nant of eigenstates of the one-particle hole Hamiltonian 



(2) 



where ho describes the host band structure and A is the 
spin splitting induced by polarized lattice ions. 

The Lowdin calculus consists in dividing the multipar- 
ticle basis states into two subsets, A and B, 



in To. The states n" are of the form M" T" , T" ^ Tq. 
To simplify the sum over n" , we approximate the diag- 
onal matrix element 'Hn"n" , which depends on both M 
and r", by the total energy of the multiparticle hole state 



r", Er 



(k,m)6r" ^Km 



It describes the interaction 



of r" with the average configuration of the ions' spins 
corresponding to the spin splitting A. We can thus write 
the Hamiltonian ^ in the following form: 
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n,n' e A . (5) 



The factor 'Hnn"'Hn"n' under the sum can be written 
as (M ® To\H\T" ® M"){M" ® T"\n\To ® M'), where 
n = M (E) To and n' = M' ® Fq. Since the denominator 
in Eq. (O is independent of M" , summing over M" is 
equivalent to inserting an identity operator, which allows 
us to write the last term as 
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(M (g) ro|?^|r")(r"|-H|ro ® W) 



Er 



Er 



We can thus treat T-L"^ as a Hamiltonian acting on ion 
states only, 

y-MM' = {M®ro\n'\To®M') 



E 

r"/ro 



{M®To\n\T"){r'\n\To®M') (6) 
Er„ — Er" 



n = 



y-AA y-AB 
y-BA y-BB 



(3) 



Set A contains all states M (E) Tq, where Fq is the P- 
particle ground state of h. Set B contains all the remain- 
ing states, in which at least one hole is excited above the 
Fermi level. We construct the effective Hamiltonian for 
the states from set A only, adding their coupling with set 
B as the second order perturbation. 
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where n,n' e A. The term ('Ho)nn' is independent 
of the ion configurations, so we set it to zero for sim- 
plicity. Thus, the effective Hamiltonian y°^ depends 
only on the ion degrees of freedom but, thanks to the 
Lowdin method, takes hole excitations into account and 
can be used to calculate the spin-wave dispersion in a 
self-consistent manner. 

The variational part of our method consists in search- 
ing for the energy E in the range where we expect to find 
the lowest eigencnergies of "H, which are the ones that we 
are interested in. For a known average spin splitting A, 
we can set E to the total energy of the hole multiparticle 



state Fo, Er„ = 



(k,m)ero 



i?k m , where i?k m is the en- 



ergy of the m-th band with wave vector k of Hamiltonian 
h, and the sum goes over all occupied eigenstates (k, m) 



Since the p-d exchange term in y, which produces the 
extra-diagonal matrix element (M ® Fo|H|F" ® M"), is 
the interaction of a single hole with an ion, the only F" 
states which have a non-zero contribution to the sum 
over F" in Eq. ([6|) are those which are created from Fq 
by just one excitation, (k, m) — >■ (k', to'). Hence, we have 
ii^ro ~ Er" = m — Eyg^i^m' and Hamiltonian ([S]) can be 
written as 



yfiM' = {M<E>To\y'\ro<E>M') 



EE 

k.k' m.m' 



f]<.,m (l /"k' ,m' ) 
£^k,m — ^k' ,m' 



(7) 



X {M E)ro\y\r"){T"\y\To ® M') 



where /k.m is the Fermi-Dirac distribution coefficient. 
The fraction in the above sum looks dangerous, as it 
may diverge in the presence of the energy bands' cross- 
ings, which would make our perturbation calculus invalid. 
However, the effective Hamiltonian for ions depends on 
the average of these factors, and will be shown immune 
to this problem. 

We write Hamiltonian ([6|) using ion spin operators and 
integrate out the hole degrees of freedom: 
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^eff ^ ^ ^ jj.g. + E E E ^ifSjsj: 

a j=l aa' j=l j' = l 



(8) 



4 



The coefficients iJJ and HJ^^ are given by 
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(9) 



wliere due to tiie condition T" ^ Tq in Eq. ([5]), for k = k' 
the summation goes over m ^ m\ and 
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A/? \ ^ \ - /k,m(l — fk',m') 



(10) 



X {s 



km.krn' '-'km'krn ^ kyji'km kmkm' / i 



where sj^^'m 



= (uk,m|s'^|Mk',m') for CT = and 
s~] = by convention. To obtain the above expres- 
sions, we substituted the hole-only part of Hamiltonian 
"Ho (H]) by the sum of P one-particle Hamiltonians /iq from 
Eq. d^l). We also used the formula (?/'k,m|^olV'k',m') = 
'5k,k'(^m,m'-E^k,m " As^^^/^,). It ariscs from the fact 
that the Bloch states 'i/'k,™ = e*'' '"Mk,m building the 
multiparticle hole states F are one-particle eigenstates 
of Hamiltonian h (U), which includes spin splitting A. 
Assuming that e*'' '' is a slowly-varying function of r 
and /(r — Rj) is constant within the unit cell around 
the j-th ion and vanishes outside of it, we obtained 
{i>Km\ri{v-n,)\4,^.^^,) = 4e*('''-'')i^^s^„kw, where 
V is the crystal volume. 

The second term of T-V^^ ([8]) describes the exchange 
interaction between the lattice ions. In the presence of 
the spin-orbit coupling it has an antisymmetric part in 
form of the Dzyaloshinskii-Moriya interaction. 
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is the antisymmetric Levi-Civita symbol 
= 1, while Mij is a pseudovector and exists 



where e^'Vo-^ 
with e 

only in the systems with broken inversion symmetry. 
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Using the Lowdin perturbation- variational calculus, 
we have thus described the problem as a lattice spin 
system coupled by exchange interaction. The low-lying 
energy states of such systems are wavelike and can be 
modeled in the small oscillations approximation. For this 
purpose, we replace the spin operators with certain non- 
linear functions of bosonic creation and annihilation op- 
erators Oj and Qj, carrying out the Holstein-PrimakofF 
bosonization4i Next, we approximate these functions 
with their power expansions around the state of satu- 
ration magnetization; 



VSa.j , 5,,- 



VSal 



^S- a]aj , (12) 



leaving in the Hamiltonian the terms which are quadratic 
in creation and annihilation operators, as only these 
terms influence the dispersion relation. This approxi- 
mation works very well in the long-wave limit, aq <^ tt, 
as the neglected magnon-magnon interactions are pro- 
portional to {aq)'^ ^ Furthermore, it justifies neglecting 
any short-range interactions between localized spins. At 
the same time, the continuous medium approximation 
becomes valid, which allows us to carry out the Fourier 
transform Cq = 7V~^/^ SjLi e^^'-^^aj. After simple alge- 
braic transformations, we arrive at the final form of the 
harmonic Hamiltonian, 
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(13) 



We call it the interaction representation as it describes 
the perturbation of the ground state by the isotropic 
Coulomb interaction (first term) and by the spin-orbit 
interaction, coupling modes of different q (remaining 
terms). The spin susceptibility of the holes is given by 
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k+q 



E\i' 



E, 



k+q,m' 



(14) 



^ '^km(k+q)m'*^(k+q)m'km ' 



where n — N/V is the density of localized spins S in the 
sample volume V and nS/S — A. The presence of the 
energy denominator shows that Xq'^ corresponds to the 
second-order part of the Hamiltonian ([7]). As promised, 
vanishing of the denominator is not harmful, due to the 
de I'Hospital rule. 

The formula (fH)) implies that Xq^ = (Xq~)* is sym- 



metric in q, while Xq i Xq^ G 7?. inherit the symmetry of 
the '0k. m eigenstates. Hence, we expect the latter to be 
symmetric with respect to q for systems which preserve 
space inversion symmetry^i^ and otherwise for systems 
which do not. The q-independent term describes the in- 
teraction of a single magnetic ion with a molecular field 
arising from the intraband spin polarization of the carri- 
ers. 
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(15) 



The corresponding term reflecting the interband polar- 



ization. 
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(16) 



arises from both Hjp ^ and the second part of Hj co- 
efficient (dni) , and cancels exactly in the full Hamiltonian 

Hamiltonian (|13p in the interaction representation de- 
scribes the spin system in terms of circularly polar- 
ized plane waves (first term), which interact with each 
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other and deform in time (remaining terms). We want 
to obtain the dispersion relation of independent, stable 
magnons. For this purpose, we diagonalize "H®* by the 
Bogoliubov transformation from aq,a]j to 6q,6q opera- 
tors, keeping in mind that we deal with the system which 
breaks the space inversion symmetry. The final form of 
the effective Hamiltonian in the spin-wave representation 
reads 

q 

where excitation modes are spin waves with dispersion 

^^^x^^^^^lWE^SFI^2, (18) 

In the case of Xq~ = X^q, fulfilled for the systems invari- 
ant under space inversion, the above formula simplifies 
to Konig's et alr^ solution, 

c.q = v'(s-xr)'-ixrp • (19) 

Furthermore, neglecting the spin-orbit coupling, when 
X^^ ~ 0, Bogoliubov transformation is unnecessary, 
and the effective Hamiltonian is already diagonalized 
by Cq, aj- operators. The negative sign of the Hamil- 
tonian (flTl) is a consequence of using the electronic con- 
vention to describe the hole-ion system. 

Apart from spin waves, which are the eigenstates of the 
stationary Hamiltonian T-L°^, dynamic excitations of dif- 
ferent physical origin may occur (e.g. Stoner spin- flips). 
They transfer a single carrier across the Fermi level to the 
state excited by the energy fko = i?k,m — ^^k-i-q,m' (see 
the denominator of Eq. (|14p . and may lead to the spin 
waves' damping. In the presence of the spin-orbit cou- 
pling, we find them likely to appear at very low energetic 
cost throughout the whole q-vector range. On the other 
hand, within the mean-field approximation and well be- 
low the Curie temperature, substitutional and thermal 
disorder are characterized by a short correlation length 
^, which ensures a well-defined spin- wave excitation spec- 
trum for q < 2tt/^. 

Figure [T] presents a typical spin- wave dispersion spec- 
trum Wq ([18]) in bulk (Ga,Mn)As calculated by the tight- 
binding computational scheme described in Sec. I. The 
hole concentration p equals 0.65 nm~'^ and the spin split- 
ting A = — 0.13eV is applied along the easy axis z fixed 
to the [001] direction by the biaxial strain Sxx = —0.()%<^ 
The latter causes a small anisotropy between Wq for spin 
waves propagating in the [001] direction and in the (001) 
plane. Furthermore, due to the lack of inversion symme- 
try in the unit-cell geometry, Xq" 7^ X-q i^i general. This 
means that, contrary to the results of the six-band k ■ p 
model^ ujq can be asymmetric with respect to the sign 
of q or A (inset). (This effect was studied experimentally 
and theoretically in the context of dielectric susceptibil- 
ity in another zinc-blende crystal, InSbJ^,) Related to 
the very small Dresselhaus spin splitting of the conduc- 
tion band, the also small ujq asymmetry is observed for 
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q VECTOR (nm"^) 

FIG. 1: (Color online) Dispersion dependence of spin waves 
propagating along [100] (equivalent to [010]) and [001] direc- 
tions in bulk (Ga,Mn)As. Inset: the Uq asymmetry due to 
the bulk inversion asymmetry of GaAs lattice modeled by the 
spds* tight-binding method, as compared to the six-band k-p 
model preserving the symmetry. Lower and upper solid lines 
denote the [412] and the opposite [4l2] propagation direction. 




q|l [010] (nm"^) -4 q || [100] (nm"^) 



FIG. 2: (Color online) Dispersion dependence of spin waves 
in two (Ga,Mn)As monolayers. The minimum shift of in the 
[110] direction to qmin = (0.43, —0.43) nm~^ can be observed. 



all but [100] and [111] q directions, and decreases with 
the growth of spin splitting. This result remains in line 
with our previous studies showing that static properties 
of (Ga,Mn)As like Curie temperature or anisotropics, to 
which ujq belongs, depend mainly on the properties of 
the carrier bands.— 

Figure [2] presents Wq in (Ga,Mn)As thin layer for p = 
0.3 nm~^ and A = —0.1 eV along the easy axis [110]. The 
host crystal consists of two unstrained infinite monolayers 
(Ga, As, Ga, As) grown in the [001] direction. The new 
effect, related to the structure asymmetry, is a shift of 
the dispersion minimum to a non-zero qmin value. 

According to the above plots, ujq for small q-vectors 
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can be described by the following general formula: 



- Wo 



(20) 



where indices fi,iy = x,y, z denote spatial directions used 
in the Einstein sum convention. The D and U constants 
are the spin- wave stiffness tensor and the Dzyaloshinskii- 
Moriya coefficient, respectively, while ujq is the spin- wave 
gap created by magnetocrystallinc anisotropics. The 
higher order terms arising from the bulk inversion asym- 
metry can be skipped as negligibly small for considered 
q- vectors. 

Alternatively, we can expand the q-dependent terms 
in LOq, considering their properties implied by Eq. p4p . 



Xc 



Ms 



where Mg ~ g^^nS is the saturation magnetization of 
the spin system. We obtain that, up to quadratic terms, 
the spin-wave stiffness in Eq. ()20|) depends only on the 
X+- term. 



D = , while U = ^tj . (22) 

Ms Ms ^ ^ 

The above constants are the subject of micromagnet- 
ics, which we shall investigate in the next section. It ig- 
nores the quantum nature of the atomic matter and uses 
classical physics in the limit of a continuous medium. 



B. Macroscopic picture 

The atomic-scale effects investigated in the previous 
section lead to the wavelike behavior of the spin system. 
In micromagnetics, these spins are replaced by classical 
vectors with their slow-varying direction n(r) described 
by the free energy functional: 
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(23) 



U''eafsn"'d^nP d^r 



where indices = x,y, z and a, /3 = x,y denote spa- 
tial and magnetization directions, respectively, and eap 
is the antisymmetric Levi-Civita symbol. The first term 
describes the anisotropy energy, which depends on the 
orientation of the magnetization with respect to the easy 
axis z. Consecutive orders of the magnetocrystallinc 
anisotropy tensor /C in principal-axis representation are 
numbered by j. The next term is the symmetric ex- 
change energy, where jA. is the exchange stiffness tensor, 
with = Ap'^. The antisymmetric part of exchange 
is expressed by the last, Dzyaloshinskii-Moriya term4— 
Since we describe the magnetization fluctuations around 



the easy axis, we do not include the derivatives of in 
the sum, as they are of higher order. 

The exchange stiffness A. can be split into two parts. 
The first one is isotropic in the magnetization direction, 
but can bear the anisotropy of the exchange interaction 
in space (e.g. caused by biaxial strain). 



(24) 



The remaining part is anisotropic with respect to the 
magnetization direction. 



(25) 



We define indices cr, (t' = +, — , referring to correspond- 
ing spin components, so that = (n^ ± irfi^j^pl and 
(21) A^^f^di^n'^d^nf^ = A^^^.d^n" d^n" . It follows from the ten- 
sors' definitions that {T^D* = Ti^':: and = 0. We 
shall use the new notation to rewrite the exchange en- 
ergy, including the Dzyaloshinskii-Moriya term, in the 
following form: 



(26) 



d^r 



In analogy to microscopic approach, we transform the 
exchange energy functional to the reciprocal space. 



T^tqi^q. a(q)ta(-q)^ + T^^g^g. a(q)«(-q) 



(27) 



We compare the above result to the effective Hamilto- 
nian (|13[) in the interaction representation. It is clear 
that A (together with U) and T correspond to coeffi- 
cients of the and Xq^ expansions in Eq. ([2T|l. respec- 
tively. Furthermore, we can identify the components of 
the above form with different physical mechanisms gov- 
erning the spin behavior. The term involving A describes 
the energy of a circularly polarized spin wave, as it arises 
from the isotropic part of exchange interaction. Hence, 
we shall call it the isotropic exchange stiffness tensor. 
The two terms involving T account for the anisotropic 
exchange, as we have chosen in Eq. ([^5]) . induced by the 
spin-orbit coupling. Hence, we shall call it the relativis- 
tic exchange stiffness tensor. Its non-zero elements imply 
that the tilting of an individual spin from the easy axis 
z to different directions has different energetic cost. As 
a consequence, the polarization of the spin wave deforms 
and acquires an elliptical shape. The linear term charac- 
terized by the constant U represents the minimum shift 
of the spin-wave dispersion dependence observed in thin 
(Ga,Mn)As layers (Fig. [2]), associated with the asymmet- 
ric exchange of Dzyaloshinskii-Moriya ([TT|) . The energy 
Wo in Eq. (j20p is related to the anisotropy constant /C in 
the full free energy functional ([^S]) . 
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We calculate the above tensors for the bulk (Ga,Mn) As 
from Fig. [1] and two monolayers from Fig. [21 For this 
purpose, we fit the coefficients of the Xq" a-^d ex- 
pansions in Eq. (pT|) on ca lnm~^ edge cube (or square) 
in q-space, centered around zero. The obtained tensors 
describe the energy of spin waves polarized in the plane 
perpendicular to the easy axis z. If they propagate in 
this plane, we call them longitudinal waves. Transverse 
spin waves propagate along z. 



1. Bulk (Ga,Mn)As 

The dispersion spectrum of spin waves propagating in 
bulk (Ga,Mn)As along two main crystal axes, [100] || x 
(equivalent to [010]) and [001] || z, is presented in Fig.[T] 
The simulated system is biaxially strained, Exx = —0.6%, 
with the hole concentration p =■ 0.65 nm~'^ and spin split- 
ting A = — O.lSeV along the easy axis z \\ [001]. 

The energy cost of exciting a circularly polarized wave 
is given by the exchange stiffness tensor 




mcV nm 



(28) 



It is expressed by a diagonal form with eigenvectors 
pointing along crystal axes, as the magnetization in the 
spins' ground state is uniform. The difference between 
its elements refiects the anisotropy of the exchange in- 
teraction in space (between the xy plane and the growth 
direction z) caused by the biaxial strain. 

In the presence of the spin-orbit coupling, the circular 
polarization can deform into an ellipse. This polarization 
anisotropy is described by the relativistic T^.+ tensor, 
which depends on the mutual orientation of magnetiza- 
tion and spin- wave propagation directions. In our system 
with the easy axis z along the spatial z direction, it takes 
the following form: 




-0.035 i 0> 
-0.020 

0; 



meV nm 



(29) 



Its zero diagonal component means that the polarization 
of transverse spin waves is circular, while the non-zero 
elements indicate that longitudinal spin waves have el- 
liptical polarization, with the shorter axis of the ellipse 
rotated to the q direction. The resulting polarizations 
are illustrated with Fig. [3]a. 

For an arbitrary propagation direction, the shape of 
the spin-wave polarization can be calculated from the 
following polarization matrix: 



/ Vxx Pxy 
\Pyx Pyy 



(30) 



(a) 



^'4 





FIG. 3: (Color online) Spin-wave polarization (exaggerated 
for clarity) , namely a shape traced out in a fixed plane by the 
spin vector rotating around the magnetization easy axis z, in 
(a) the bulk (Ga,Mn)As from Fig. [1] with z \\ [001] and (b) 
after changing z to [010]. 



where Pa/s = A^pq^q,y, or in more detail 



Pii =q(A + 9^eT++)qT 
q(A-5HeT++)qT 



Pyy 



-xt 



Pxy ^Pyx^- 3mqT++q'^ - 3m Xo'' ■ 
The ellipse is the solution of the equation 
n°'n^Pafj = const , 



S, (31) 



(32) 



where and are the in-plane components of the 
magnetization unit vector n. The ratio of its longer and 
shorter main axes, a and h, can be derived from the eigen- 
values of the polarization matrix pop as 



a/b 



Pxx + Pyy + \J iPix - PyyY + ^Ply 
\ Pxx + Pyy - \J iPxx ~ PyyY + ^Ply 



(33) 



Kinetically, a longitudinal elliptical spin wave can be 
viewed as a circular one, which experiences the Lorentz 
contraction in the direction of motion, h = a^Jl — w^/c^, 
traveling with the velocity v, where c is the speed of 
light. Figure |4] presents this spin-wave relativistic veloc- 
ity in the relevant range of q vectors along [100] and [110] 
directions. 

As already mentioned, T is induced by the spin-orbit 
coupling, which connects the symmetries of the lat- 
tice with spins. Thus, similarly to magnetocrystalline 
anisotropics, it depends on the magnetization direction 
with respect to the crystal axes. The T tensor calcu- 
lated in the analyzed (Ga,Mn)As system with the easy 
axis z changed to the [010] direction would acquire the 
following form: 




meV nm 



(34) 
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0.08 




q VECTOR (nm ) 

FIG. 4: (Color online) Relativistic velocity of spin waves prop- 
agating in two crystallographic directions: [110] (dotted line) 
and [100] (solid line), as a function of the q- vector magnitude. 

which reveals its dependence on the biaxial strain. (An 
appropriate rearrangement of the matrix elements gives 
the T++ tensor for z along [100].) Now the longitudi- 
nal spin waves propagate in the xy \\ (010) plane and 
experience the 'Lorentz contraction' as described above, 
while the transverse waves propagate along the new easy 
axis z. Their polarizations are additionally deformed 
by the potential of the strained crystal. The latter is 
stretched equally in the [100] and [010] directions and 
compressed in the [001] direction, and so are the polar- 
izations. The magnitude of their deformations is given 
by the q-independent terms of the polarization matrix 
p (jSO]) . Figure [3]b is illustrative of this effect. 

Technically, the relativistic exchange might lead to 
macroscopic anisotropics in the system, if they did not 
average out for spin waves propagating in different di- 
rections. However, one can imagine a weak anisotropy 
arising from the described phenomena in asymmetrically 
shaped samples, where the largest number of similarly 
polarized modes can exist along the longest dimension. 
It would then be a candidate for an explanation of the 
weak uniaxial anisotropy of the [100] and [010] crystal 
axes observed in some (Ga,Mn)As samples4^ 

2. (Ga,Mn)As layers 

In thin layers of (Ga,Mn)As (Fig. [2]), we observe the 
minimum shift of the spin-wave dispersion to a non- 
zero Qmin valuc, which was not present in bulk. It is 
a hallmark of the Dzyaloshinskii-Moriya asymmetric ex- 
change (jlip . arising from structure inversion symmetry 
breaking. The mechanism of this interaction is demon- 
strated in our numerical simulations for the two mono- 
layers of (Ga,Mn)As from Fig. [2j 

The structure of the modeled lattice is shown in Fig. [5] 
Thanks to the lack of inversion symmetry at the mid- 
point of each Ga-As bond, the exchange interaction of 
Dzyaloshinskii-Moriya is allowed, Ujj/ • Sj x Sj' . Accord- 




FIG. 5: (Color online) Two infinite monolayers of GaAs. The 
Ga and As sites are numbered by their elevation in the growth 
direction [001]: I's are situated at (•, -,0), 2 at (■, ■, jo), 3 at 
(•, ■, I a) and 4 at (•, •, |a). The molecular field of magnetic 
Mn ions is introduced according to Sec. [H The strongest in- 
teractions between the nearest neighbors are marked by solid 
bonds. The created solid lines define mirror planes in the zinc- 
blende crystal, (110) and (110). Since the analyzed structure 
has two-dimensional periodicity, the carrier momenta in the 
growth direction are quantized and spin waves propagate in- 
plane only. 



ing to Moriya's rulcs;^ the u^- vector of each atom pair, 
indicated by an arrow at each bond, is perpendicular to 
their mirror plane. Its sense is always the same when 
we go from Ga to As atom. Since in our theoretical ap- 
proach outlined in Sec. Ill Al wc have used the Holstein- 
Primakoff transformation to describe spin waves as small 
fluctuations around the magnetization direction z, only 
the terms contribute to the spin-wave dispersion re- 
lation. For the above reasons, we expect the maximum 
effect of the Dzyaloshinskii-Moriya interaction in systems 
magnetized along the u vectors: perpendicular to [110] 
or [110]. First, we consider the magnetization z set along 
the [110] direction and a spin wave propagating perpen- 
dicular to it, q II [110]. The dotted line is the wave- 
front, along which all spins must be in phase. Along this 
wavefront, the constant u vector (big arrows) tilts the 
spins perpendicular to itself and to each other to mini- 
mize the energy of the Dzyaloshinskii-Moriya interaction. 
The Ujji vector is a function of the distance between the 
spins, Rj — Rj' . Hence, when we move from one spin pair 
to another in a regular lattice structure, their chirality 
Sj X Sj' minimizing the energy is constant in magnitude 
and antiparallel to Ujy . Therefore, each spin will be ro- 
tated with respect to its neighbors by a constant angle 
around a constant axis. The rotation direction is parallel 
or antiparallel to the spin-wave propagation direction q, 
depending on whether xijji is parallel or antiparallel to z. 
In this way, a cycloidal structure with a period A forms 
in the spin system. Since the modulation occurs along 
the [110] direction, this is where we observe the disper- 
sion minimum shift by g„iin = 27r/A (Fig. [2). For the 
magnetization z pointing along the [110] direction, the 
u vectors (small arrows) cancel out when looking along 
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this direction, and no frustrated structure of lower en- 
ergy can arise. Since the Dzyaloshinskii-Moriya interac- 
tion operates in the sample plane, neither will we observe 
its hallmarks when z is perpendicular to this plane. 

For the quantitative analysis of the described effects, 
we fit the spin- wave dispersion presented in Fig. [5] with 
the form ((20|) . In Fig. [6] (full circles) we sweep the mag- 
netization z in the sample plane (001) and report the ob- 
tained angle between z and qmin together with the energy 
gain Wo — Wq^^,;^ due to the Dzyaloshinskii-Moriya inter- 
action, and the magnitude of qmin- From simple algebra 
we have q,nin = -^D^^U and wq - Wq„,i„ = q^f^jj^ D q^in. 
As expected, the energy gain is the largest when z and 
qmin are perpendicular to each other, for z || [110], and 
drops to zero for z || [110] (Fig. [6]a,b). The first ar- 
rangement is accompanied by the strongest frustration 
of the spin system, qmin = 0.6 nm^^ (Fig. Elc), resulting 
in the spin cycloid with period A = 10 nm. In the other 
arrangement, no spin frustration arises and the collinear 
spin order is preserved. This behavior accounts for an in- 
plane anisotropy with the easy and hard axes along the 
[110] and [110] directions, respectively. It is easy to no- 
tice that if we mirror-reflect the sample (or cquivalently, 
move the bottom Gal layer to the top Ga5 position), the 
two axes switch — [110] will become the hard axis and 
[110] the easy axis (Fig. [6l open circles). Consequently, 
the sign of qmin and the orientation of the spin cycloid 
will change. The described anisotropy arises from the 
net Dzyaloshinskii-Moriya interaction. Although it is a 
surface effect, it should not be confused with the surface 
anisotropy, which docs not occur in the zinc-blende (001) 
surface. 

In the above reasoning, it is important to remember 
that we deal with the hole-mediated ferromagnetism. 
While the magnetic lattice ions substitute only Ga sites 
and would seem oblivious to the inversion asymmetry of 
the Ga-As pairs, the carriers interact with all surround- 
ing atoms. Thus, the system of magnetic ions interacting 
through the hole carriers, as described by Eq. ([8]), is sen- 
sitive to all symmetry properties of the crystal. 

As we have shown above, the Dzyaloshinskii-Moriya 
interaction in thin (Ga,Mn)As layers leads to the frus- 
tration of spins in the ground state. While in bulk 
(Ga,Mn)As the diagonal form of tensor A (pS)) depicts 
the uniform ground state magnetization, the exchange 
stiffness tensors in thin layers reveal a spin cycloid ori- 
ented along one of the mirror planes: 



[110] 



[110] 



[1T0] 



A 
T 



0.93 
-0.28 



-0.28 
0.93 



mcV nm 



0.09- O.Olt 
0.05 + 0.005i 



0.05 
0.09 



0.005i 
0.006i 



(35) 



meV nm 



Their eigenvectors point along the highest symmetry 
axes, [110] and [110]. This effect is equivalent to ap- 
plying strain along one of these directions, which would 
account for a magnetoelastic anisotropy, like the one in 
biaxially strained samples. However, contrary to the 







O" 


7C/2 


-V- 









> 




0) 




E, 


0.4 


E 




3^ 


0.2 


1 

o 




3 











0.6 


T 
E 
c 


0.4 


"e 


0.2 


cr 







' (a) 















Till 7C 3/27t 

MAGNETIZATION DIRECTION z 



27C 



FIG. 6: (Color online) Parameters of the minimum shift qmin 
in two (Ga,Mn)As monolayers from Fig. [2] as a function of 
the magnetization direction z for the analyzed system (full 
circles) and for its mirror-reflection (open circles), a) The 
direction of the qmin vector with respect to z, b) energy gain 
[Jo — '^q„,i„5 c) magnitude of the minimum shift qmin. 



case of a uniformly strained sample, all effects related to 
the Dzyaloshinskii-Moriya interaction vanish in the bulk 
limit. 

The A, T and U constants let us think of the mag- 
netic system in terms of the interactions by which it 
is governed, as described by the Hamiltonian in 
the interaction representation. Alternatively, if we want 
to deal with independent stationary magnons (with al- 
ready deformed polarization), represented by the Hamil- 
tonian (|17p . we calculate the spin- wave stiffness D, re- 
lated to A by the formula ([22]). The difference between 
these two pictures is especially apparent at finite tem- 
peratures, where the length of the magnetization vector, 
Af(r), decreases and the formula becomes 



D(r) = ^A(T) 
^ ' M{T) ^ ' 



(36) 



The higher the temperature, the stronger we have to tilt 
the magnetization vector in order to excite the quantum 
of spin waves (a magnon). Hence, A(T) decreases with 
temperature faster than D(T). 



C. Spin-wave stiffness and Curie temperature 

This section discusses the relation between the spin- 
wave stiffness and Curie temperature, based on which 
we define the normalized spin- wave stiffness parameter, 
Dnor- Depending only weakly on the hole density p and 
spin splitting A, it makes a convenient tool for experi- 
mentalists to estimate the spin-wave stiffness values D 
given the Curie temperature and the magnetization of 
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the sample. We provide quantitative numerical results 
on Unor for bulk (Ga,Mn)As obtained in the spds* tight- 
binding computational scheme (see Sec. I), and compare 
them to the outcome of the previously employed six-band 
k ■ p model We also clarify that the surprisingly large 
magnitude of D in these systems results from the p-like 
character of the periodic part of Bloch function. 

The standard relation between spin-wave stiffness D 
and Curie temperature Tq in a cubic crystal, which is 
often used in the literature on the topic, reads 



D = 



2 

nn 



2(3 +1) 



(37) 



where r^n is the nearest neighbor distanccj^ Derived for 
short-range interactions, it is interesting to find out how 
it is modified when considering the actual nature of the 
spin-spin exchange. This question has a number of im- 
plications. For instance, both carrier relaxation time, 
which is limited by magnon scattering, and the quantita- 
tive accuracy of the mean-field approximation grow with 
D/Tq, as the density of spin waves at given temperature 
T diminishes when D increases. It is worth noting that 
a simple adaption of this formula to (Ga,Mn)As by re- 
placing Tun with {xnQ)~^/'^ resulted in an overestimation 
of the p-d exchange integral by an order of magnitude 

Allowing for spatially modulated structures, the mag- 
netic ordering temperature Tq for a carrier-controlled 
ferromagnet is given by a solution of the mean-field 
equatioufii^ 



/3'x(q,7^)X5(q,T) = l 



(38) 



where x ^nd xs are the carrier and lattice ion spin sus- 
ceptibilities. In the simple case of a parabolic band they 
are proportional to the Pauli and Curie- Weiss magnetic 
susceptibilities, respectively. 

First, we consider the two-band model of carriers re- 
siding in a simple parabolic band, where all that is left 
of the A. tensor of Eq. (|23|) is the scalar isotropic ex- 
change stiffness A. According to the previous section, 
it is related to the quadratic coefficient of the expan- 
sion x(g) « Xo + nSP'^Cq^, hy D = nSP^C = 2A/nS. 
Additionally, we assume that the ground state of the sys- 
tem corresponds to the uniform ferromagnetic ordering, 
q = 0, and take xs{(l,T) in the Curie form. 



Using Eqs. 



nS{S + l) 
^^^^'^)= 3fcBT 

and ([5^ . we obtain 

SfceTc CiT -> 0) 



D = 



(5 + i)x(o,r = rc) 



(39) 



(40) 



If the values of both spin splitting A at T — s- and /cbTc 
are much smaller than the magnitude of the Fermi energy 
\E-f\, the carrier susceptibility is given by 



x(.)4p(i^p)F(^ 



(41) 



The total density of states at 
p{E-p) = m* kp / {nh)'^ , where m* 
mass, and 



F 



Ef = h^k^l{2m*) is 
is the carrier effective 



^2kpJ 


1 


k-p 




2ky + q 




^ 2 ~ 






2kp — q 



= 1 - 



(42) 



is the Lindhard function. We obtain from these equations 

kBTc 



D 



4(5- 



(43) 



We see that, in agreement with the notion that magnetic 
stiffness increases with the range of the spin-spin interac- 
tion, Tnn in Eq. ^ is replaced by l/ikpV^) in Eq. ((43)) . 
which scales with the range of the carrier-mediated in- 
teractions. Indeed, according to the Ruderman-Kittel- 
Kasuya-Yosida theory;^ the magnitude of the ferromag- 
netic exchange integral decays at small spin-spin dis- 
tances r as l/{rkp) and reaches the first zero at r sa 
2.2/fcF. 

Since in semiconductors l/fcp ^ aq, the above for- 
mulae imply that D/Tq is rather large in systems with 
carrier-controlled ferromagnetism. A question arises as 
to how the ratio D/Tq would be modified, if the com- 
plexities of the valence band were taken into account. 

As already noticed by Konig et al.^ the values of ex- 
change stiffness for (Ga,Mn)As are greater in the six- 
band model with a spin-orbit coupling than in the case 
of a simple parabolic band. As argued by these authors, 
due to the multiplicity of the valence bands, the Fermi 
level for a given carrier concentration is much lower than 
in the two-band model, and hence both carrier polariza- 
tion and exchange stiffness are greaterJ^ On the other 
hand, Brey and Gomez-Santos^i assign large values of D 
to higher magnitudes of Tq in the multi-band model. 

To check these suggestions we plot in Fig. [7] the values 
of dimensionless parameter 



4(5"- 



l)klD 



kuTi 



(44) 



for various hole concentrations p and fcp = (37r^p)^/'^. 
The ratio D/Tq for bulk (Ga,Mn)As is computed us- 
ing the spds* tight-binding (Sec. |l| and the previously 
employediS six-band k ■ p model of the semiconductor 
band structure. Guided by the results of the two band 
model, one could expect -Dnor ^ 1 in the limit of small 
spin polarizations, A <C I^'fI- In contrast to these ex- 
pectations, Dnor is of the order of eleven and, moreover, 
weakly dependent on the polarization, altered here by 
changing the hole concentration and the spin-orbit split- 
ting. Furthermore, the experimentally observed biaxial 
strain magnitudes have only slight effect on the stiffness 
tensor (see Fig. [1] and Eq. ([28]) .— This indicates that the 
single parameter Dnor can serve to estimate the magni- 
tudes of D and A if only the Mn magnetization AI and 
hole concentration p are known. 
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FIG. 7: (Color online) Normalized spin-wave stiffness Dnor 
in bulk (Ga,Mn)As as a function of the hole concentration 
p and the spin-orbit coupling (inset; the value of the spin- 
orbit coupling in (Ga,Mn)As is indicated by dotted line) for 
different spin splittings A. Filled and empty markers indicate 
the results of the spds* tight-binding and the six-band k ■ p 
model, respectively. 



Knowing that neither the spin-orbit coupling nor the 
multiplicity of carrier bands can explain the large spin- 
wave stiffness, we turn our attention to the matrix el- 
ements (Mk.m|s'^|Mk+q.?Ti') in thc Spin susceptibility of 
holes (fTl)) . Neglecting the spin-orbit coupling, Uk,m is 
a product of the spin part Sm and the real spatial part 
W'k. m- Thus, we can write 



{U]i,m\S |Mk+q,i 



,m I'^k+q,^ 



for (T = +, — . In the parabolic two-band model with 
its carrier wave functions described by plane waves 
'0k,m(i") = ^^e"*^"", the periodic part Uk,™ = 
Hence {sm.\sm') = Smm' and iuk,m = and F(q) = 
x(q)/x(0) simplifies to the Lindhard function (|42|l . More 
realistic models take into account the periodic lattice po- 
tential, which mixes different atomic orbitals and leads to 
the k-dependent modulation of Uk.m- (The eventual com- 
position of hole states in the spds* tight-binding model 
is presented in Fig. [5]a-c.) From this it follows that 
the scalar products ('ii'k,m|'Wk+q,m') are q-dependent, and 
can be only smaller than in the simple model. As a conse- 
quence, F(q) has a steeper slope, as presented in Fig.[8]d. 
This explains why the magnitude of the stiffness tensor 
is so large in (Ga,Mn)As and related ferromagnets. Ac- 
tually, the fact that the p-type character of the carrier 
wave functions affects in this way the g-dependence of 
dielectric susceptibility has been already noted4^i^ 
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FIG. 8: (Color online) Atomic composition of heavy, light 
and spin-orbit split-off hole bands (a, b and c panels) as a 
function of the fc- vector. It is given by the squares of norms of 
projections of the carrier state on atomic s, p, d or s* orbitals 
of As or Ga atoms (solid and dotted lines). Panel d presents 
the Lindhard function F{q) calculated in the tight-binding 
model (p = 0.65 nm~^ and A = — 0.13eV, with the spin-orbit 
coupling set to zero) and in the parabolic band model with 
the same Fermi energy = — Q.27eV. Dashed lines are their 
square form fits. 



D. Spin waves' contribution to magnetization 

In this section, we address the problem of applicability 
of the mean- field model to the description of temperature 
dependence of magnetization in the analyzed systems. 
While it is a known fact that the critical fiuctuations do 
not change much the mean-field Curie temperature Tc in 
the presence of long-range exchange interactions,— sev- 
eral papers^"— discuss how it is lowered by spin waves. 
Also, the influence of dilution and disorder on spin waves 
and the Curie temperature, which can be crucial in sam- 
ples with low Mn content, was studied by many authors 
with varying conclusions.—^— In response, we propose a 
simple model which takes into account the correct num- 
ber of spin-wave excitations and thermal depolarization 
of spins. However, we argue that spin waves also do not 
lower the mean-field Tc due to thermal decoherence of 
quantum spin system. 

The mean-field approximation employed in our model 
allows for a reasonable overall description of ferromag- 
netism in the analyzed systems. It reduces the problem 
of lattice spins coupled by the exchange interaction to 
that of noninteracting spins in the molecular field S (jlSp . 
Their magnetization is described by the self-consistent 
equation. 



M{T) = A/(0) 



^s(A(r),r) 

fcRT 



(45) 



The Brillouin function B5 ignores the actual nature of 
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thermal fluctuations and their correlations, and assumes 
that every spin fluctuates independently, which in our 
theory corresponds to the high-g limit of spin-wave exci- 
tations (see Fig. [T|). At low temperatures, however, the 
long- wavelength magnons of much lower energies can ex- 
ist. Since each of them carries a moment of /ie, they yield 
a strong contribution to the temperature dependence of 
magnetization 



5.58 



Af(r) = il/(0)(l-(iV5)-i^(nq)T 



(46) 



where {nq)T is the thermal average of spin- wave excita- 
tions in each mode. The latter can be modeled by the 



Bose-Einstein distribution. 



MIT 



l/(eXp(L^q/fcBr) 



1), as spin waves are bosons. Then, replacing the sum- 
mation in Eq. (|46p by an integral and putting ojq « Dq^, 



one obtains the well known T'^/^ Bloch lav 



dq (nq)T = C3/2 



3/2 



kBT 
D 



3/2 



(47) 



where C3/2 ~ 2.612 is the Riemann zeta function. 

Remembering about the spin-wave gap created by 
magnetic anisotropy, which allows for ferromagnetism 
in low dimensional systems despite the Mermin- Wagner 
K loq + Dq^. (We neglect the 



theoremi^ we have 



Dzyaloshinskii-Moriya coefficient, which vanishes in bulk 
limit.) In that case, the above law is modified to^ 



dq 



MIT 



= L 



13/2 



-u,Q/kuT\ 3/2 



kBT 
D 



3/2 



(48) 



where Li3/2 {e~'^°/^'^'^^ is dc Jonquicrc's function. 

Furthermore, the Bose-Einstein statistics allows for the 
unlimited number of spin waves in each mode (on the 
other hand, Konig et ali^ assume that there can be only 
2S spin waves of each q, which is too strict). In fact, their 
total number cannot be larger than 2NS, corresponding 
to complete magnetization reversal. We handle this by 
introducing a fictions mode with zero energy, which is 'oc- 
cupied' by the spin waves which have not been excited in 
reality. In this way we can treat the problem with classic 
Bose-Einstein condensate methods)^ the total number 
of bosons occupying all modes is always 2NS, their zero- 
energy mode constitutes the 'condensate' phase, while 
the excited spin waves constitute the 'thermal cloud'. 
The total number of spin waves in the limit of infinite 
crystal volume is therefore given by 



^q/f^nT 



vq<qD 



Additionally, while at zero temperature the system is 
described by a pure state 4'o, where exciting a spin wave 
costs the energy 
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FIG. 9: (Color online) Temperature dependence of magne- 
tization for bulk (Ga,Mn)As from Fig. [U calculated by the 
mean- field model (gS]), T^/^ Bloch law ^ and its modified 
version H48|) . and the 'spin- wave condensate' (|51JK The exper- 
imental dependence for the sample from Ref. (2^ is analyzed 
in Sec. [ml 



at non-zero temperatures it is described by a mixture of 
pure states with a certain number of spins flipped 
by one-particle thermal excitations, p ~ X]nP".l^")(^"l- 
Thus, its magnetization drops to M(T) according to the 
mean- field Brillouin function (|45|) . The energy cost of 
exciting a spin wave is now given by 



^p„[(6^vi/„|H°^|6tvI/„)-(*„|H'=«|*„)] . (50) 



Since the spin-wave dispersion (fT8|) depends on tem- 
perature T almost exclusively via the spin splitting A and 
is approximately proportional to it, we can estimate the 
above expression by Wq "Pr. 



'q L^n t-n M{0) ~ A/(0) ~ "^q A(0) ' 

where M„ is the magnetization in the state V^n. Spin 
waves are thus perturbations of the thermal state of lat- 
tice spins and not of the ground state, which we can 
model just by replacing M(0) with M(T) (or A with 
A(T)) in the dispersion relation Wq. They additionally 
lower the magnetization to M'{T). Remaining in the 
limit of small oscillation approximation, we obtain the 
following set of equations: 



M(T) 



A(T) 



M{T) 
M'{T) 



: M(0)Bs 

: M{T) 



'SE{A'iT)y 



— min 




(51) 



,2M{T) 



where 



j(A(r), T)/kBT. The spin splitting A(T) 



is induced by the mean field of the lattice spins described 
by the first equation, while A'(r) additionally takes into 
account their depolarization due to spin waves. 
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Figure [9] presents the temperature dependence of mag- 
netization as described by the above methods for the bulk 
(Ga,Mn)As from Fig. □ The T^/^ Bloch law and its 
modified version apply to low temperatures only. We see 
that the values of D determined neglecting the anisotropy 
gap would be overestimated. The 'spin-wave conden- 
sate' includes the magnon contribution to magnetization 
in the whole temperature range, and leads to signifi- 
cantly lower Curie temperature than in the mean-field 
model. This method uses the correct bound for the num- 
ber of spin-wave excitations. It also naively attempts to 
solve the problem of the well-known shortcomings of the 
spin- wave theory^ introduced by the Holstein-Primakoff 
transformation (fT2|l . as it includes thermal disorder by 
depolarizing the lattice spins with temperature. How- 
ever, quantum-mechanical intuition suggests that the 
spin waves should vanish completely at higher temper- 
atures. Non-zero temperature leads to the loss of infor- 
mation about the system: pure states are replaced by 
mixed states, namely the thermal states e~^^/Tre~^^, 
and the quantum correlation of spins vanishes. This 
can be reflected in the p-d Hamiltonian by replacing 

■ s, with /3E,j(4)(^;|) = which is 

simply the mean-field model. Similarly, in the vicinity of 
the Curie transition, where the problem of critical fluc- 
tuations arises, the mean-field models of systems with 
long-range exchange interactions are known to work very 
wefl4^ 

The next section provides further insight into the ap- 
plicability of the presented models by comparing their 
results to experiment. 



III. COMPARISON TO EXPERIMENTAL DATA 

As already mentioned in Sec. I, the theory of ex- 
change stiffness developed within the six-band p-d Zener 
modeU^ describes quantitatively the width of stripe do- 
mains in (Ga,Mn)AsJ^ Recently, Gourdon et al^ car- 
ried out a detailed analysis of the magnetic domain struc- 
ture and magnetic properties of an annealed 50 nm-thick 
Gao.93Mno.07As layer with a perpendicular magnetic easy 
axis and the Curie temperature of 130 K. Two employed 
experimental methods yielded an upper and lower limit 
of the isotropic exchange stiffness A as a function of tem- 
perature T. As found by these authors from examining 
the domain- wall velocity, the higher values of A{T), de- 
termined from the lamellar domain width, are reliable. 

We model the sample in the tight-binding computa- 
tional scheme for bulk (Ga,Mn)As (see Sec.|T|. In order to 
determine the material parameters for numerical simula- 
tions, we start by estimating the effective Mn content Xcs 
from the measured low-temperature spontaneous magne- 
tization M{T — > 0) = 39kA/m. Taking into account the 



hole contribution, Mc 



5kA/m,— implies the mag- 



netization of the Mn spins Ms = 44kA/m. This value 
corresponds to the effective Mn content XeS = 4.3% and 
the spin splitting A = — 0.13eV. No direct measurements 
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FIG. 10; (Color online) Theoretical reconstruction of experi- 
mental data on exchange stiffness j4 as a function of temper- 
ature T obtained in Ref. i22| from the analysis of the lamellar 
domain width in Gao.93Mno.07As. 



of the hole concentration are available for this sample, so 
we estimate its magnitude from the effective and total 
Mn content, XcS and x = 7%. Assuming that intersti- 
tial Mn donors had survived the annealing process and 
formed antiferromagnetic pairs with the substitutional 
Mn acceptors^^^ we obtain p ~ (3a;off/2 — x/2)no = 
0.65 nm-3,^ 

From Fig. [71 for the given value of p and A we 
find Dnor ~ 10.5, which gives the spin- wave stiffness 
D = l.lmeVnm^ {A = 0.21 pj/m) at T -)> OK. Know- 
ing this, we can calculate the temperature dependence 
of magnetization according to Sec. IIIDI (Fig. [9|). In the 
mean- held picture (solid lines), the magnetization of lat- 
tice spins Ms is described by the Brillouin function (|45p . 
The magnitude of hole magnetization \ Mc\ decreases with 
temperature proportionally to Ms, according to Eq. (jlSp . 
The resulting magnetization M{T) is compared to the ex- 
perimental curve (circles). We obtain a good agreement 
with the measured data, especially near the Curie tran- 
sition, and Tc = 127 K. At low temperatures, we plot 
the outcomes of the T^/^ Bloch law (H7|l and its modified 
version (|48p employing the calculated spin- wave stiffness 
value D (dotted lines). The modified Bloch law, adjusted 
to include the spin-wave gap, gives very good agreement 
with the experimental dependence, which indicates that 
in this regime the spin- wave excitations are solely respon- 
sible for demagnetization. Near the Curie transition, we 
reconstruct the measured Curie temperature and mag- 
netization values with the Brillouin function, which sug- 
gests that the temperature destroys the spin- wave coher- 
ence and recalls the mean-field picture. 

To reconstruct the A{T) trend obtained from the 
magnitudes of lamellar domain width,— we again make 



use of the fact that the susceptibility Xq~ dl) de- 
pends on temperature almost exclusively via spin split- 
ting. Thus, and according to the formula the ex- 
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FIG. 11: (Color online) Theoretical reconstruction of experi- 
mental data on exchange constant J for samples with varying 
effective Mn content ajeji, as determined from the tempera- 
ture dependence of magnetization using the T^l'^ Bloch law 
(large crosses) and from the Curie temperature (small crosses) 
in Ref. [l7|. Empty circles indicate the region where the ex- 
perimental Tc exceeds the theoretical predictions. Inset: ex- 
perimental temperature dependence of magnetization for the 
sample with Xeff = 3.3% fit by the standard Bloch law with 
J = 0.26 meV (Ref. [13), the modified Bloch law (0 with the 
theoretically obtained J = 0.12 meV and by the mean-field 
Brillouin function. 



change stiffness scales with temperature as A(r)^ We 
use the experimentally determined M(T) and the cal- 
culated A{T — > 0) = 0.21 pj/m to estimate the ex- 
change stiffness values for the remaining temperatures 
as A{T) = A{0)M{T)/M{0). As shown in Fig. [lOl this 
procedure correctly reproduces the experimental A{T) 
trend. However, the same value of A{T — >■ 0), which 
has been shown to successfully describe the experimental 
M{T) dependence with the modified Bloch law, is twice 
as large as the exchange stiffness constant determined 
from the measured lamellar domain width. 

Potashnik et ali^ evaluated the isotropic exchange 
constant J for the set of optimally annealed (Ga,Mn)As 
layers with varying Mn content. The values derived from 
the temperature dependence of magnetization, using the 
standard Bloch law, and from the Curie temperature 
within the three-dimensional Heisenberg model?^ were 
similar. 

Our theoretical reconstruction of the experiment by 
the spds* tight-binding model of bulk (Ga,Mn)As is 
demonstrated in Fig. 1111 To obtain the presented results, 
we have estimated the effective Mn content Xcs from the 
measured low-temperature magnetization M[T 0), 
and then increased it by about 10% to include the hole 
contributioui^ Assuming that annealing removed all in- 
terstitials and each remaining substitutional Mn pro- 
duces one hole carrier, the hole density p = Xeff n-o- For 
the obtained values of p and XeS we have calculated 
the mean-field Curie temperatures Tq and the exchange 



constant J = D{xcf[no)'^^^ /2S assuming that Mn ions 
form a cubic lattice. We note that this form is equiv- 
alent up to a few percent to that employed in Ref. [13: 
J = i:)(47rx7io/24)2/3/2S', where X is the total Mn con- 
tent (N. Samarth, private communication). Since for the 
samples with low Mn content our values of Tc are much 
lower than the experimental ones, we conclude that the 
corresponding estimates of J (indicated by empty cir- 
cles) are not reliable. On the other hand, we reconstruct 
Curie temperatures for samples with Xcs > 2.5% (filled 
circles) , but this time the theoretical values of J are much 
smaller than the experimental ones. This discrepancy 
points to the importance of the anisotropy-induced en- 
ergy gap in the spin-wave spectrum. As illustrated in 
Fig. [11] (inset) for the sample with XcS = 3.3%, fitting 
the experimental M{T) curves with the standard Bloch 
law neglecting the gap leads to higher values of J than 
those expected theoretically. At the same time, the mod- 
ified Bloch law (jlS]) employing the theoretical exchange 
constant J = 0.12 meVnm^ reconstructs the analyzed 
M{T) trend. We notice that it perfectly describes the 
mild slope of the low-temperature M(T), contrary to the 
standard Bloch law Similarly to the case of the Gour- 
don et al. experiment;^ the mean-field model works very 
well at higher temperatures. 

In a series of experiments, the spin- wave stiffness was 
determined by examining spin precession modes excited 
by optical pulses^^ and under ferromagnetic resonance 
conditions ji^Ti^i According to these works, the experi- 
mental findings are strongly affected by gradients of mag- 
netic anisotropy, presumably associated with carrier de- 
pletion at the surface and interface, which also affect the 
character of spin pinning. We also note that no infiuence 
of the magnetic field on the hole spins, visible as a de- 
viation of the Landc factor from the value g = 2^ was 
taken into account in the employed Landau-Lifshitz equa- 
tions. With these reservations we show in Fig.[T2|the ex- 
perimentally evaluated values of D plotted as a function 
of the nominal Mn concentration x for various as-grown 
and annealed samples of Gai_a;Mnj;As. These findings 
are compared to the results of ab initio computations^^ 
(dashed line) and our theory for the hole concentration 
p = xno and p = 0.3a:no (solid lines). When comparing 
theoretical and experimental results, one should take into 
account that the actual Mn concentration Xcn is smaller 
than X, particularly in as-grown samples. As seen, our 
theory describes properly the order of magnitude of the 
spin-wave stiffness D but cannot account for a rather 
large dispersion in the experimental data. 

In thin (Ga,Mn)As layers, described in Sec. Ill B 21 
the Dzyaloshinskii-Moriya interaction can be observed in 
form of a spin- wave dispersion minimum shift. It leads 
to the formation of a cycloidal spin structure and uniax- 
ial in-plane anisotropy of the [110] and [110] directions, 
with the easy axis determined by the sample geometry. 
Both these phenomena have been a subject of consider- 
able interest in recent years<^"— While the long-period 
spin structures have not been hitherto observed, the uni- 
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FIG. 12: (Color online) Compilation of theoretical results 
(lines), obtained with no adjustable parameters, and exper- 
imental data (markers) on the spin-wave stiffness D as a 
function of the nominal Mn concentration x in Gai_2,Mna; As. 
Solid symbols denote annealed as well as etched and hydro- 
genated samples fRef. [2ll) . Empty symbols denote as-grown 
samples. The value g = 2 was used to convert D in magnetic 
units (Tnm^) to energy units (meVnm^). 



axial anisotropy is commonly present in (Ga,Mn)As.— 
However, it was shown by gradual etching^ or by inves- 
tigating different samples^ to be insensitive to the layer 
thickness. This is not the case for the anisotropy caused 
by the Dzyaloshinskii-Moriya interaction, which is a sur- 
face effect and vanishes with increasing layer thickness. 
To test our theory, one should pick very thin high-quality 
samples grown in the [001] direction, which additionally 
excludes the anisotropy of the surface, with the easy axis 
along one of the in-plane diagonal directions. The cy- 
cloidal spin structure could be then observed e.g. under 
a magnetic force microscope or by neutron scattering. 
The sample should be probed along the [110] and [110] 
directions to find the long-period modulation of magneti- 
zation only along the one perpendicular to the easy axis. 



ties have been expressed by the phcnomenological pa- 
rameters of micromagnetic theory. We have noticed that 
the strength of ferromagnetic order described by the 
isotropic exchange stiffness is significantly amplified by 
the p-like character of carrier wave functions, as com- 
pared to the simple parabolic band model. Furthermore, 
we have found various effects reflecting the tendency of 
the spin-orbit interaction to pervade every aspect of car- 
rier dynamics. They produce the relativistic corrections 
to spin- waves given by the anisotropic exchange stiffness 
tensor and the asymmetric Dzyaloshinskii-Moriya coeffi- 
cient. The latter accounts for the cycloidal spin arrange- 
ment and the accompanying uniaxial in-plane anisotropy 
of diagonal ([110] /[1 10]) directions in thin layers, which 
can result in a surface-like anisotropy in thicker films. 
Quantitative results on the stiffness constant have been 
provided in form of a normalized parameter, which as- 
sumes the value Dnor ~ 1 1 over a wide range of Mn and 
hole concentrations in (Ga,Mn)As. They agree with the 
previous k ■ p calculationsiSiii but predict significantly 
smaller values of spin- wave stiffness than those resulting 
from ab initio computations^^ 

The above theories have been applied to analyze the 
related experimental data on the stiffness parameter and 
the temperature dependence of magnetization. Our ba- 
sic theoretical model has not managed to reconstruct all 
stiffness values obtained by various experimental meth- 
ods. In all cases we have reconstructed the entire range of 
magnetization dependence on temperature. At low tem- 
peratures, it can be understood within the the modified 
Bloch law^ employing the values of spin-wave stiffness 
calculated by our model. At higher temperatures, the 
mean-field theory becomes justifiable owing to thermal 
decoherence and the long range character of spin-spin in- 
teractions. At the same time, we have reproduced only 
partly the stiffness values obtained from analyzing pre- 
cession modes in (Ga,Mn)As thin films. Our results may 
allow to separate bulk and surface effects, as well as bring 
to light the pining phenomena and the role of inhomo- 
gencities in experiments examining precession modes in 
slabs of carrier-controlled ferromagnetic semiconductors. 



IV. SUMMARY 

We have investigated spin waves and exchange stiff- 
ness in thin layers and bulk crystals of ferromagnetic 
(Ga,Mn)As described by the spds* tight-binding com- 
putational scheme. Using the proposed variational- 
perturbation calculus, we have described the analyzed 
systems and their spin-wave excitations. Their proper- 
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